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ABSTRACT 


On  2  April  1997,  strong  winds  blew  through  the  central  coast  of  California 
that  were  accompanied  by  an  intense  jet  streak  and  upper  front.  The  event  was 
analyzed  with  standard  synoptic-scale  DIFAX  charts  and  mesoscale  charts  for 
comparison.  The  mesoscale  model  used  was  the  Coupled  Ocean/ Atmosphere 
Mesoscale  Prediction  System  (COAMPS)  developed  by  Naval  Research  Laboratory 
(NRL),  Monterey,  California. 

CO  AMPS  captured  features  that  were  not  diagnosed  on  the  synoptic  charts. 
Height,  isotach  and  temperature  fields  showed  more  detail,  although  observations 
were  not  sufficient  to  completely  verily  the  model’s  level  of  detail.  CO  AMPS  was 
heavily  influenced  by  the  topographic  field  modeling  lee  troughs  and  mountain 
waves  along  the  Sierra  Nevada  mountains.  A  strong  mountain  wave,  initiated  by 
the  upper  front,  ocurred  in  central  California  during  this  time  period.  The  wave 
troughs  correlated  to  wind  maxima  at  the  surface,  including  one  near  San 
Francisco  Bay  where  winds  as  high  as  66  mph  were  reported. 
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I.  INTRODUCTION 


The  headlines  in  three  San  Francisco  Bay  area  newspapers  on  the  morning  of 
3  April  1997  all  pertained  to  strong  winds  that  blew  through  the  region  on  the  2nd. 
The  San  Francisco  Examiner  read  "Windblown  Opener",  describing  the  San  Francisco 
Giants  opening  day  game  on  the  2nd.  The  Monterey  Herald's  "High  winds  blow  out 
the  lights"  discussed  the  loss  of  power  to  more  than  6000  Monterey  county  homes 
caused  by  the  strong  winds.  And  the  San  Jose  Mercury  News  reported  winds  as 
strong  as  66  mph  on  Mount  Diablo  (east  of  Oakland  in  the  Diablo  mountain  range) 
under  the  headline  "The  Answer,  My  Friend...".  These  reports  accompanied  with  the 
knowledge  of  the  general  synoptic  situation  at  that  time  which  depicted  a  strong  jet 
streak  in  the  area  (see  Chapter  IV  for  details),  provided  the  elements  for  follow-on 
work  proposed  by  Pauley  et  al.  (1996)  regarding  jet  streaks,  upper  fronts  and  their 
influence  on  surface  winds.  This  case  study  examines  this  strong  low-level  wind 
event. 

The  objective  of  the  study  is  to  compare  the  Coupled  Ocean/Atmosphere 
Mesoscale  Prediction  System  (COAMPS)  analyses  of  the  upper  front/jet  streak 
structure  with  available  observations  and  with  standard  DIFAX  charts  and  to  diagnose 
the  behavior  of  the  analyzed  upper  front/jet  streak  and  its  effect  on  the  low-level 
winds.  The  study  will  also  explore  mountain  waves  depicted  by  COAMPS  that  were 
generated  as  the  upper  front  interacted  with  the  Sierra  Nevada  range,  east  of  the  Bay 
area. 
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A.  UPPER  FRONTS 

Upper  fronts  are  formed  by  secondary  circulations  about  jet  streaks  as  a  result 
of  geostrophic  confluence  and  cross-frontal  geostrophic  wind  shear  as  described  in 
Shapiro  and  Keyser  (1990),  and  Keyser  and  Shapiro  (1986).  Figure  1  is  a  schematic 
depiction  of  these  dynamic  processes,  which  provide  the  framework  for  analyzing  the 
secondary  direct  and  indirect  circulations  responsible  for  jet  streak  propagation  and 
upper  frontogenesis.  The  secondary  circulations  caused  by  geostrophic 
confluence/d  iff  I  uence  without  thermal  considerations  (Fig.  la)  are  indirect  in  the  exit 
region  of  the  jet  streak  and  direct  in  the  entrance  region.  The  indirect  circulation 
forces  warm  air  to  sink  acting  to  increase  the  frontal  thermal  gradient  and  vorticity 
while  the  direct  circulation  forces  warm  air  to  rise  acting  to  decrease  the  frontal 
thermal  gradient  and  vorticity.  Horizontal  shear  modifies  the  location  of  the 
secondary  circulations  when  thermal  advection  is  present  and  also  acts  frontolytically 
on  the  cyclonic  shear  side  of  the  jet  streak  (Fig.  1b)  and  frontogenetically  on  the 
anticyclonic  shear  side  when  there  is  an  along-jet  positive  thermal  gradient.  A 
combination  of  the  two  is  more  typical  with  cold  advection  likely  west  of  the  trough 
axis,  as  shown  schematically  in  Fig  2c.  Figures  2d-f  from  Keyser  and  Shapiro  (1986) 
also  shows  three  other  orientations  and  their  resultant  shift  of  the  induced  secondary 
circulations. 

There  have  been  several  models  of  upper  front  structure  during  the  past  60 
years  (the  reader  is  directed  to  Keyser  and  Shapiro  (1986)  for  a  complete  historical 
and  dynamical  review).  Figure  3  from  Keyser  and  Shapiro  (1986)  is  a  diagram  of  the 
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current  model  for  analyzing  upper  fronts  and  tropopauses.  The  cyclonic  shear  zone 
of  the  upper  front  is  on  the  order  of  100  km,  much  less  spread  out  than  earlier  theory 
suggested.  This  structure  allows  for  the  exchange  of  air  parcels  between  the 
troposphere  and  stratosphere  and  also  possibly  the  downward  momentum  transport 
of  kinetic  energy  to  low  levels  where  boundary  layer  processes  can  mix  the  energy 
to  the  surface  (Pauley  et  al.  1996).  Tracers  such  as  ozone  (O3)  and  potential  vorticity 
(PV)  are  used  to  locate  and  identify  upper  fronts.  Because  PV  is  an  easily  derived 
quantity  from  current  numerical  weather  prediction  models,  it  is  normally  used  for 
upper  front  analysis  and  will  be  used  in  this  case  study.  Figure  4  from  Pauley  et  al. 
(1996)  is  a  cross-section  of  an  upper  front  and  jet  streak  with  PV  contours  delineating 
the  upper  front. 

B.  MOUNTAIN  WAVES 

When  stably  stratified  air  is  forced  to  rise  over  a  topographic  barrier,  a 
disturbance  is  created  that  is  carried  away  from  the  source  region  by  gravity  waves 
(Durran  1990).  By  convention,  when  mountains  are  the  topographic  barrier 
involved,  the  waves  are  referred  to  as  mountain  waves.  High  amplitude  mountain 
waves  can  also  be  associated  with  strong  surface  winds  that  blow  down  along  the 
lee  slope  of  the  mountain  (Durran  1990)  i.e.,  downslope  winds.  Key  elements  for 
this  phenomena  to  occur  are  that  the  wind  direction  be  within  30"  of  perpendicular 
to  the  mountain  ridge  line  and  that  the  upstream  temperature  profile  exhibits  a  layer 
of  strong  stability  near  mountaintop  level  (Durran  1990). 
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A  case  study  of  a  January  1972  windstorm  in  Boulder,  Colorado  was  also 
presented  by  Durran  (1990).  Figure  5  shows  the  isentropic  cross-sections  for  this 
event.  In  Fig.  5a,  the  atmospheric  conditions  are  modeled  using  actual  upstream 
(negative  values  on  the  x-axis)  conditions  that  contain  a  strong  stable  layer  near 
mountaintop  level.  This  condition  generates  a  strong  mountain  wave  signature 
downstream  from  the  mountain.  In  Fig.  5b,  the  stable  layer  is  removed  from  the 
model  run  with  all  other  parameters  the  same.  Notice  the  absence  of  the  wave 
signature  downstream  from  the  mountain  indicating  that  the  stable  layer  interaction 
with  the  mountain  initiated  the  wave.  In  the  case  presented  in  this  thesis,  the  upper 
front  provided  the  stable  layer  for  this  interaction. 

The  thesis  is  organized  as  follows.  The  COAMPS  data  assimilation  system 
used  in  this  research  is  briefly  described  in  Chapter  II,  to  the  extent  necessary  for 
basic  understanding.  The  methodology  presented  in  Chapter  III  identifies  the 
FORTRAN  program  used  to  run  and  plot  the  data  fields  from  the  model,  followed  by 
a  synoptic  overview  in  Chapter  IV  usirig  standard  charts  to  describe  the  broad  scale 
atmospheric  conditions  during  the  case  study.  A  comparison  of  COAMPS  model  data 
to  the  standard  chart  data  is  then  conducted  in  Chapter  V  followed  by  a  discussion 
in  Chapter  VI  on  the  gravity/mountain  waves  depicted  by  COAMPS  in  this  case. 
Finally  results  of  the  study  and  recommendations  for  follow-on  work  are  presented 
in  Chapter  VII. 
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II.  COAMPS  OVERVIEW 


COAMPS,  the  Coupled  Ocean/Atmosphere  Mesoscale  Prediction  System, 
consists  of  an  atmospheric  nonhydrostatic  model  and  a  hydrostatic  ocean  model. 
This  section  will  provide  a  brief  description  of  the  atmospheric  model,  the  only 
component  of  COAMPS  used  in  this  study.  The  COAMPS  atmospheric  model  is 
designed  for  either  idealized  or  real-data  simulations.  This  case  study  uses  real-data 
assimilation  which  consists  of  data  quality  control,  analysis,  initialization  and  forecast 
components  which  are  briefly  described  below.  For  a  detailed  account  of  the  model, 
the  reader  is  directed  to  Hodur  (1997),  the  source  for  this  overview. 

A.  GRID  CONFIGURATION 

Polar  stereographic,  Mercator,  spherical  grid  or  Lambert  conformal  projections 
are  available  with  the  model.  An  Arakawa  C  grid  on  a  Lambert  conformal  projection 
is  used  for  this  case  study.  The  model  utilizes  nested  grids,  with  up  to  nine  meshes 
allowed  provided  there  is  a  3:1  reduction  in  grid  spacing  between  grids.  For  this 
case  study  a  two-mesh  nested  grid  was  used,  81  km  for  the  coarse  grid  and  27  km 
for  the  fine  grid.  Detailed  analysis  was  conducted  on  the  27  km  inner  grid.  The 
model  uses  the  sigma-z  vertical  coordinate  scheme  consisting  of  30  levels  with 
approximately  9-10  levels  in  the  planetary  boundary  layer. 

B.  INITIAL  CONDITIONS 

Several  surface  fields  must  be  initialized  before  running  COAMPS.  Terrain 
height  is  obtained  from  the  U.S.  Navy  20'  resolution  terrain  field  (approx.  20  km)  or 
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the  Defense  Mapping  Agency  (DMA)  Digital  Terrain  Elevation  Data  (DIED)  level  1 
data  (100  m  resolution).  The  U.S.  Navy  20'  terrain  field  was  used  for  this  case  study 
as  seen  in  Fig.  6.  Surface  albedo,  roughness  and  ground  wetness  are  also  specified 
primarily  by  using  monthly  climatological  global  fields  bilinearly  interpolated  to  the 
COAMPS  grids.  Ground  temperature  is  obtained  from  a  surface  energy  budget  over 
land  and  NOGAPS  analyzed  surface  temperature  (held  constant  for  model  integration 
cycle)  over  water. 

Extensive  quality  control  (QC)  programs  (Baker  1992)  are  used  to  check 
observational  data  for  errors.  The  data  that  pass  the  QC  checks  are  then  mapped  to 
the  model  grids  using  a  multivariate  optimum  interpolation  (MVOI)  analysis 
technique  (Lorenc  1 986  and  Barker  1 992).  While  processing  data  for  this  case  study, 
it  was  discovered  that  there  was  a  problem  with  the  COAMPS  implementation  of  the 
MVOI  scheme  for  the  inner  mesh.  Briefly,  the  MVOI  scheme  was  performed  using 
four  volumes  for  the  inner  mesh  that  did  not  overlap.  First-order  discontinuities  at 
the  edge  of  the  analysis  volumes  were  seen  in  this  case,  because  the  western  two 
volumes  used  the  data-sparse  algorithm,  while  the  northeastern  volume  used  the  data- 
rich  algorithm.  When  the  volumes  were  assembled,  derivative  quantities  such  as 
divergence  had  large  values  at  the  interface  between  analysis  volumes,  bisecting  the 
domain  shown  in  Fig.  6  in  both  x  and  y  directions.  To  solve  this  problem,  the 
volumes  were  made  to  overlap  so  that  each  region  sampled  had  a  similar  data 
distribution  resulting  in  a  more  continuous  analysis  increment  field  (Holt,  personal 
communication).  Observations  used  in  this  study  were  obtained  from  radiosondes. 
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pibals,  AIREPS,  ACARS,  SSM/1,  surface  obs,  cloud  track  winds,  DMSP  and  NOAA 
satellites. 

An  incremental  or  full  update  cycle  can  be  used  for  the  analysis.  In  the 
incremental  update  cycle,  which  was  used  for  this  study,  the  analyzed  increment 
fields  are  interpolated  from  the  16  mandatory  pressure  levels  from  1000  to  TO  mb  to 
the  model's  sigma-z  vertical  levels  and  added  to  the  model's  sigma-z  level  first  guess, 
preserving  the  model's  vertical  structure.  Once  the  initial  fields  are  obtained,  they 
must  be  initialized  using  a  variational  method  to  reduce  the  generation  of  spurious 
high-frequency  oscillations  and  to  assure  that  vertical  propagating  sound  waves  are 
not  generated. 

C.  FORECAST  MODEL 

Lateral  boundary  conditions  for  the  COAMPS  coarse  grid  are  obtained  from 
the  most  recent  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS) 
forecast.  The  NOGAPS  fields  are  horizontally  interpolated  using  a  bicubic  spline  to 
the  COAMPS  coarse  grid,  then  interpolated  in  the  vertical  to  the  model  vertical 
coordinate.  The  grids  utilize  one-way  interaction,  i.e.,  outer  to  inner  for  interior 
boundary  conditions.  Second-order  numerical  differencing  is  used  throughout  the 
model  except  for  diffusion  which  is  fourth-order.  The  time  differencing  scheme  used 
is  split  explicit. 

Clouds  and  precipitation  are  treated  differently  depending  on  the  grid  spacing. 
Below  10  km  resolution,  COAMPS  has  explicit  moist  physics  (Rutledge  and  Hobbs 
1983)  and  performs  as  a  cloud  model.  Above  10  km,  as  in  this  case,  COAMPS  uses 
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explicit  treatment  of  the  moist  physics  for  nonconvective  clouds  and  precipitation  and 
the  Kain  -  Fritsch  (Kain  1990,  Kain-Fritsch  1993)  cumulus  parameterization  for 
convective  clouds  and  precipitation.  The  Louis  et  al.  (1982)  scheme  is  used  to  model 
surface  fluxes,  the  Harshvardan  et  al.  (1987)  method  is  used  for  atmospheric  radiation 
and,  the  Deardorff  (1980)  parameterizations  for  subgrid-scale  mixing  is  used. 

D.  COAMPS  SPECIFICS  FOR  APRIL  1997  CASE 

For  the  case  presented  in  this  thesis,  COAMPS  was  run  in  data  assimilation 
mode  using  a  6-h  incremental  update  cycle,  beginning  0000  UTC  1  April  1997.  A 
NOGAPS  analysis  was  used  for  the  background  at  the  initial  time,  after  which  a  6-h 
COAMPS  forecast  was  used  as  background.  A  two-nest  version  was  used  with  61  x 
61  grid  points  on  the  outer  81  km  resolution  grid  and  73  x  73  grid  points  on  the 
inner  27  km  resolution  grid,  both  centered  at  35°N  120°W  on  a  lambert  conformal 
projection.  Thirty  vertical  levels  were  used. 
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III.  METHODOLOGY 


The  sigma-level  COAMPS  analyses  were  post-processed  to  pressure  levels  at 
a  50  mb  interval  from  1050  mb  to  100  mb  and  displayed  using  FORTRAN  code  that 
was  developed  by  Prof.  P.  M.  Pauley  (Naval  Postgraduate  School)  and  Dr.  E.  H. 
Barker  (Naval  Research  Laboratory,  Monterey)  for  previous  case  studies  including  the 
1991  dust  storm.  The  code  was  revised  to  accept  COAMPS  data  rather  than 
NORAPS  (Navy  Operational  Regional  Atmospheric  Prediction  System)  data.  Namely 
the  vertical  coordinate  definition,  the  grid  dimensions,  and  the  input  routines  needed 
to  be  revised  to  use  the  COAMPS  fields.  The  plotting  portion  of  the  program  is 
based  on  National  Center  for  Atmospheric  Research  (NCAR)  graphics. 

This  diagnostics  package  can  compute  sea  level  pressure  using  several 
different  algorithms.  The  Mesinger  (1990)  technique  was  chosen  as  being  less  noisy 
than  the  other  possibilities.  In  this  technique,  temperatures  are  interpolated 
horizontally  in  below-ground  regions  by  setting  the  Laplacian  of  the  temperature  field 
to  zero  at  below-ground  points  and  using  the  neighboring  above-ground  points  on 
a  given  pressure  level  as  boundary  conditions.  The  solution  was  obtained  by 
sequential  overrelaxation.  The  resulting  temperatures  were  used  to  compute  heights 
or  pressure  surfaces  at  below  ground  points.  The  sea-level  pressure  values  were  then 
obtained  by  extrapolating  downward  from  the  nearest  pressure  surface  that  is  above 
sea-level. 
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The  data  for  each  date-time-group  were  transferred  using  ftp  from  Naval 
Research  Laboratory  (NRL)  in  Monterey,  California  and  then  unpacked.  The 
diagnostics  code  was  run  separately  for  each  date-time-group  with  several  choices 
available  for  output.  Plan-view  charts  with  conventional  fields  such  as  heights, 
winds,  vertical  motions,  and  temperatures  were  primarily  used.  Vertical  cross- 
sections  were  also  plotted  and  used  for  diagnosing  theta  surfaces,  potential  vorticity, 
mountain  waves,  and  the  interaction  of  the  upper  front  with  topography. 
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IV.  SYNOPTIC  OVERVIEW 


This  chapter  describes  the  synoptic  situation  that  encompasses  the  high  wind 
event  in  central  California  on  2  April  1997.  National  Center  for  Environmental 
Prediction  (NCEP)  ETA  DIFAX  charts  at  standard  levels  will  be  used  to  describe  the 
synoptic-scale  features'  initial  state  and  evolution  during  the  period. 

A.  0000  UTC  1  APRIL 

A  high  amplitude  longwave  pattern  exists  from  850-300  mb  over  the  eastern 
North  Pacific  through  the  Northern  Plains  of  the  United  States  at  this  time  (Figs.  7- 
10).  The  ridge  axis  of  interest  to  this  case  study  is  centered  at  approximately  140  W 
extending  well  into  the  Gulf  of  Alaska.  The  trough  axis  of  interest  lies  along  120  W 
over  eastern  Washington  and  Oregon  with  the  base  of  the  trough  over  central 
California  and  southern  Nevada.  A  second  ridge  is  centered  at  95  W  along  the 
western  Minnesota  border. 

At  300  mb  (Fig.  7),  observations  indicate  a  65  m/s  jet  streak  is  slightly 
downstream  of  the  ridge  axis  while  45  m/s  winds  extend  southward  to  the  northern 
California  border.  Note  that  the  eta  analysis  is  too  weak,  with  a  jet  streak  maxima 
of  only  110  kt,  approximately  55  m/s.  Although  a  closed  low  has  not  yet  formed  in 
the  base  of  the  trough,  a  -35°  C  cold  pocket  is  located  over  central  and  eastern 
Oregon  at  500  mb  (Fig.  8)  with  weak  cold  air  advection  into  the  trough  base. 
Moderate  warm  air  advection  is  occurring  at  the  ridge  axis. 
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At  the  surface  (Fig.  1 1),  a  1035  mb  high  over  the  eastern  north  Pacific  and  a 
1004  mb  low  located  in  southern  Nevada  is  causing  a  moderate  pressure  gradient 
over  California  with  the  strongest  gradient  in  the  southwestern  part  of  the  state.  A 
wavy  cold  front,  supported  by  the  shortwave  at  850-700  mb  (Figs.  9  and  10),  extends 
southwest  from  northeastern  Montana  through  southern  California. 

B.  1200  UTC  1  APRIL 

During  the  following  12  h,  the  ridge  builds  and  moves  east  about  5° 
developing  a  slight  tilt  to  the  northeast.  The  trough,  however,  is  quasi-stationary  and 
does  not  deepen.  The  65  m/s  jet  streak  at  300  mb  (Fig.  12),  verified  by  the  Salem, 
Oregon  rawinsonde  observation,  has  propagated  down  the  coast  to  a  location  over 
western  Washington  and  British  Columbia  aligned  north  to  south.  The  -35°  C  cold 
pocket  at  500  mb  (Fig.  13)  has  also  moved  slightly  to  the  east  over  western  Idaho. 
700  mb  shows  no  significant  features  at  this  time  or  for  the  remaining  valid  times  and 
is  therefore  omitted  for  the  rest  of  the  overview. 

At  850  mb  (Fig.  14)  a  lobe  of  the  ridge,  oriented  west  to  east,  has  built 
eastward  through  the  central  Rocky  Mountain  region  causing  the  height  contours  to 
lie  west  to  east  and  generating  flow  normal  to  the  Sierra  Nevada  range.  The  trough 
associated  with  the  embedded  shortwave  remains  through  the  southwest  along  with 
a  quasi-stationary  isotherm  pattern,  signaling  a  stalled  front  in  the  region.  Height  falls 
at  San  Diego  are  20  m  at  this  level  but  only  1 0  m  at  500  mb  (Fig.  1 3),  reflecting  the 
low-level  frontal  feature. 
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The  surface  high  is  quasi-stationary  building  to  1038  mb  while  the  low  filled 
to  1008  mb  and  moved  into  northwestern  Arizona  (Fig,  15).  The  isobars  become 
oriented  northeast  to  southwest  over  California  with  an  inverted  trough  analyzed  by 
NCEP  over  eastern  California  and  are  more  evenly  spaced  than  12  h  earlier.  The 
wavy  front  is  quasi-stationary  in  the  southwest  but  pushed  farther  east  in  the  central 
U.S,  by  the  west  to  east  ridge  lobe  described  at  850  mb. 

C.  0000  UTC  2  APRIL 

The  ridge  continues  to  build  to  the  northeast  and  the  trough  deepens  slightly, 
gaining  a  much  more  prominent  northeast  to  southwest  tilt  during  this  12  h  period. 
The  second  ridge  builds  sharply  back  to  the  northwest,  enhancing  the  blocking 
pattern.  The  65  m/s  jet  streak  at  300  mb  (Fig.  16)  centered  over  central  Oregon  also 
becomes  oriented  northeast  to  southwest  with  50  m/s  winds  extending  to  Oakland. 
Observations  to  confirm  the  strength  of  the  jet  streak  were  not  available. 

At  500  mb  (Fig.  1 7)  a  5440  m  closed  low  is  analyzed  over  southern  Nevada, 
coincident  with  the  -35°  C  cold  pocket.  The  temperature  gradient  becomes  stronger 
over  central  California  while  heights  begin  to  fall  significantly  in  southern  California 
with  San  Diego  reporting  a  70  m  fall  during  the  past  12  h. 

The  shortwave  trough  through  the  southwest  at  850  mb  (Fig.  18)  deepens 
slightly  but  remains  quasi-stationary  along  with  the  isotherm  pattern.  Enhanced  lee 
troughing  is  evident  through  central  California.  Height  falls  range  from  19  m  at  San 
Diego  to  44  m  at  Flagstaff. 
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The  1038  mb  high  at  the  surface  (Fig.  19)  moves  northeast  and  the  low  over 
northern  Arizona  along  the  quasi-stationary  wavy  front  deepens  to  1003  mb, 
producing  a  somewhat  tighter  pressure  gradient.  Some  troughing  is  evident  along  the 
San  Joaquin  Valley,  a  reflection  of  the  lee  side  trough  at  850  mb  (Fig.  18). 

D.  1200  UTC  2  APRIL 

The  longwave  pattern  has  changed  little  in  the  previous  12  h.  The  trough 
moves  slightly  to  the  southwest  and  is  now  anchored  by  the  newly  formed  8920  m 
closed  low  at  300  mb  (Fig.  20)  over  extreme  southwestern  Nevada.  The  strongest 
reported  wind  at  300  mb  is  57  m/s  out  of  the  north-northeast  at  Oakland.  The  jet 
streak  is  analyzed  accordingly  for  speed  and  is  depicted  as  having  elongated  and 
begins  to  wrap  around  the  closed  low. 

At  500  mb  (Fig.  21),  180  m  height  falls  occur  at  San  Diego  as  the  closed  low 
moves  southwestward,  becoming  nearly  vertically  stacked  with  the  closed  low  at  300 
mb  and  lower  levels. 

Deep  height  falls  are  also  taking  place  at  850  mb  (Fig.  22)  over  the  southwest. 
San  Diego  shows  a  41  m  height  fall  for  the  past  12  h  as  the  closed  low  becomes 
entrenched  in  the  region  and  is  coincident  with  the  upper  levels.  The  lee  side  trough 
through  central  California  is  enhanced  as  the  gradient  tightens  from  north  to  south. 

The  quasi-stationary  surface  high  has  weakened  to  1037  mb  but  the  quasi- 
stationary  low  along  the  wavy  front  has  deepened  to  1000  mb,  thus  maintaining  the 
relative  isobaric  packing  (Fig.  23).  The  inverted  trough  appears  to  be  closer  to  the 
coast,  rather  than  over  the  Central  Valley  area. 
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E.  0000  UTC  3  APRIL 

A  significant  change  in  the  longwave  pattern  can  be  detected  at  this  time.  The 
second  ridge  to  the  east  begins  to  break  down  as  the  previously  quasi-stationary 
shortwave  gains  energy,  moves  to  the  northeast  and  suppresses  the  second  ridge. 
The  300  mb  pattern  around  the  closed  low  remains  virtually  unchanged  as  the  low 
slowly  drifts  to  the  southwest  (Fig.  24).  A  55  m/s  jet  streak,  aligned  as  before,  is 
analyzed  by  NCEP,  supported  by  a  52  m/s  report  from  Oakland.  A  weaker  45  m/s 
jet  streak  is  depicted  southeast  of  the  closed  low,  possibly  related  to  the  jet  streak 
propagating  around  the  low. 

The  closed  low  at  500  mb  (Fig.  25)  also  drifts  to  the  southwest  and  fills  by  20 
m.  Height  rises  are  found  throughout  the  entire  region.  Although  Oakland  is  still 
reporting  strong  45  m/s  northeasterly  winds,  the  tightest  temperature  and  height 
gradient  is  located  well  off  the  southern  California  coast,  in  an  area  of  no  upper  air 
observations. 

The  isotherms  associated  with  the  quasi-stationary  shortwave  at  850  mb  have 
moved  through  the  area  to  the  southeast  (Fig.  26).  They  have  been  isolated  from  the 
shortwave  by  a  bubble  high  that  has  split  away  from  the  eastward  protruding  lobe 
of  the  offshore  ridge  that  was  discussed  earlier.  The  height  contours  have  realigned 
slightly  to  the  northeast,  weakening  the  lee  side  trough. 

The  surface  pressure  gradient  has  eased  in  southern  California  as  the  low 
along  the  wavy  front  has  moved  to  the  northeast  over  northwest  New  Mexico  and 
filled  to  1004  mb  (Fig.  27).  A  weak  inverted  trough  is  analyzed  from  southern  to 
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central  California  along  the  valley  region.  The  1037  mb  high  remains  in  place  and 
is  keeping  the  gradient  tight  in  northern  California. 
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V.  COAMPS  COMPARISON 


Although  the  operational  NCEP  charts  such  as  shown  in  the  previous  chapter 
are  produced  from  48  km  eta  model,  the  fields  are  highly  smoothed  before  producing 
the  DIFAX  graphics.  Consequently,  the  27  km  mesh  of  the  COAMPS  model  has 
more  detail  than  the  standard  NCEP  charts.  Height  fields  and  isobars  show  more 
small  amplitude  features  but  were  in  general  agreement  as  to  location  of  the 
contours.  Wind  fields  show  more  mesoscale  structure  than  are  typically  analyzed  at 
synoptic  scales,  and  other  parameters  such  as  temperature  also  show  similar  levels 
of  detail  (not  shown).  Higher  resolution  observations  were  normally  not  sufficiently 
available  to  verify  this  level  of  mesoscale  detail,  leaving  questions  as  to  what  was  real 
and  what  features  were  model  artifacts.  The  biggest  difference  between  the  synoptic 
and  mesoscale  analyses  was  the  topographic  influence  for  COAMPS,  namely, 
COAMPS  depicted  mountain  waves  and  small-scale  lee  troughs  that  were  not 
apparent  on  the  NCEP  charts. 

A.  TIME  SERIES  COMPARISON 

The  300  mb  and  surface  charts  were  chosen  for  this  comparison  because  they 
are  the  essence  of  this  study.  300  mb  provides  the  best  level  to  examine  the  jet 
streak  and  the  surface  represents  the  meteorological  effects  that  are  best  recognized. 

1.  0000  UTC  1  April 

Figure  28  shows  the  300  mb  COAMPS  analysis  chart  with  heights  as  solid 
lines,  isotachs  are  shaded  in  increments  of  5  m/s  starting  at  30  m/s,  and  observations 
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are  also  in  m/s  with  flags  equal  to  25  m/s,  a  full  barb  equal  to  5  m/s,  and  a  half  barb 
equal  to  2.5  m/s  (Note:  this  convention  is  maintained  throughout  the  paper.). 
Comparing  this  chart  to  the  NCEP  chart  (Fig.  7),  one  can  see  the  obvious  differences 
in  the  representation  of  the  isotachs.  NCEP  shows  a  smooth,  elongated  90  kt 
(approximately  45  m/s)  jet  streak  extending  well  into  central  California  whereas 
COAMPS  depicts  a  more  detailed  pattern  with  the  45  m/s  isotach  barely  reaching  the 
northern  California  coast.  Rawinsonde  observations  from  Oakland  (37.5  m/s)  and 
Salem  (37.5  m/s)  seem  to  verify  COAMPS  depiction. 

The  COAMPS  surface  chart  (Fig.  29)  shows  good  agreement  with  NCEP  (Fig. 
11)  regarding  isobar  location  and  general  surface  pressure  pattern.  Isobars  are  solid 
lines,  isotachs  are  shaded  as  before  but  begin  at  15  m/s  vice  30  m/s,  and 
observations  are  in  m/s  as  before.  There  is  one  noteworthy  difference  however. 
COAMPS  shows  a  small-scale  trough  forming  in  the  lee  of  the  Sierra  Nevada  range 
in  north-central  California,  while  NCEP  shows  straight  isobars  with  no  troughing 
evident. 

2.  1200  UTC  1  April 

The  irregular  isotach  pattern  is  again  evident  on  the  COAMPS  300  mb  chart 
(Fig.  30).  The  maximum  isotach  of  55  m/s  (approximately  1 10  kt)  is  just  touching  the 
northern  California  border  whereas  the  NCEP  110  kt  isotach  extends  almost  as  far 
south  as  San  Francisco  (Fig.  12).  Because  the  Salem  rawinsonde  is  missing  wind 
direction  and  speed,  verification  of  the  COAMPS  maximum  isotach  is  not  possible. 
It  also  appears  that  the  30  m/s  isotach  extending  south  into  the  Bay  area  on  the 
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COAMPS  chart  forced  the  110  kt  isotach  farther  south  than  it  should  have  on  the 
NCEP  chart  due  to  smoothing. 

The  surface  charts  are  beginning  to  show  disagreement  in  the  isobar  pattern 
over  land.  Note  the  small-scale  troughs  that  COAMPS  is  depicting  in  the  lee  of  the 
Sierra  Nevada  range  and  along  the  coast  (Fig.  31).  NCEP  shows  broad  inverted 
troughing  extending  along  the  coast  with  a  manually  analyzed  "trof"  (Fig.  15)  in  the 
approximate  location  of  the  COAMPS  lee  trough.  This  suggests  that  COAMPS  is 
successfully  analyzing  a  feature  that  is  known  to  be  present  but  not  properly 
represented  by  the  synoptic-scale  DIFAX  chart. 

3.  0000  UTC  2  April 

The  models  are  generally  in  agreement  as  to  the  location  and  intensity  of  the 
jet  streak  at  this  time.  Both  COAMPS  (Fig.  32)  and  NCEP  (Fig.  16)  show  the  55  m/s 
(approximately  110  kt)  isotach  over  north-central  California  with  stronger  winds 
embedded  upstream.  Although  NCEP  depicts  a  slightly  stronger  jet  maximum  (130 
kt)  than  COAMPS  (60  m/s)  ,  observations  are  not  sufficient  or  available  to  verify 
either.  Of  particular  interest  however,  is  the  onset  of  a  small-scale  trough  that  is 
analyzed  by  COAMPS  in  the  lee  of  the  Sierra  Nevada  range.  The  trough  is 
accompanied  by  a  distinct  wavy  pattern  of  the  isotach  field  in  the  same  region 
suggesting  possible  correlation  of  these  two  features. 

The  surface  differences  are  the  same  as  1200  UTC  1  April  except  that  NCEP 
(Fig.  19)  is  no  longer  indicating  a  "trof"  in  the  lee  of  the  Sierra  Nevada  range.  On 
the  COAMPS  chart  (Fig.  33),  three  separate  wind  maxima  are  being  depicted  in 
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northern  California  with  two  over  the  land  and  one  broad  maximum  off-shore.  The 
two  maxima  over  land  seem  to  be  located  in  conjunction  with  the  two  small-scale 
troughs  that  were  identified  earlier.  The  next  valid  time,  1200  UTC  2  April,  proves 
to  be  the  most  interesting  and  dynamic  in  the  sequence  and  will  be  the  primary  focus 
for  the  remainder  of  this  paper. 

B.  UPPER-LEVELS  AT  1200  UTC  2  APRIL 

Starting  at  300  mb,  compare  Fig.  20  the  NCEP  chart  with  Fig,  34,  the 
COAMPS  chart  for  the  same  valid  time,  1200  UTC  2  April.  The  NCEP  chart  shows 
a  smooth,  elongated  1 10  kt  jet  streak  extending  through  central  California.  On  the 
COAMPS  chart,  note  again  the  more  detailed  pattern  in  the  isotach  field  and  the  two 
distinct  jet  maxima  of  65  m/s,  approximately  10  m/s  stronger  than  NCEP  has 
analyzed. 

However,  the  observational  network  is  not  dense  enough  to  verify  the  strength 
or  the  location  of  these  jet  maxima.  Mesoscale  details  such  as  these  come  from  the 
six-hour  COAMPS  forecast  that  was  used  as  background  for  the  analysis.  The  Reno 
rawinsonde  observation  is  within  the  isotach  maximum  at  the  California/Nevada 
border,  but  has  no  wind  observation.  The  assumption  is  that  the  windspeed  is  strong 
and  the  balloon  was  blown  away.  Also  note  that  sometimes  even  where  there  are 
observations,  significant  differences  between  the  analysis  and  the  observations  exist. 
The  Oakland  rawinsonde  observation  has  a  windspeed  of  57.5  m/s  but  is  inside  the 
60  m/s  isotach.  The  opposite  situation  occurs  at  Vandenburg  which  has  an 
observation  of  60  m/s  but  is  just  inside  the  50  m/s  isotach.  However,  the  latter 
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observation  is  approximately  100  km  west  of  aircraft  observations  of  27.5  and  17.5 
m/s  in  a  location  where  the  analysis  is  between  30  and  35  m/s.  Unable  to  determine 
if  these  observations  are  used  in  the  analysis,  one  can  only  conclude  that  the  analysis 
possibly  has  problems  depicting  this  large  horizontal  shear.  The  NCEP  analysis 
portrays  values  of  approximately  53  m/s  and  38  m/s,  respectively,  at  these  two 
locations  which  yields  an  even  weaker  horizontal  shear  than  COAMPS.  In  regions 
with  weaker  wind  shear,  the  COAMPS  analysis  draws  more  closely  to  the 
observations,  e.g.,  Medford,  Oregon. 

At  500  mb,  for  the  same  valid  time,  NCEP  has  analyzed  straight  contours 
through  north-central  California  and  Baja,  California  (Fig.  21).  Compare  this  with  the 
COAMPS  chart  (Fig.  35),  which  shows  small-scale  troughing  of  the  height  contours 
in  both  north-central  and  Baja,  California,  regions  of  considerable  topographic 
variability  (Fig.  6).  At  this  level,  the  rawinsonde  observations  at  Oakland, 
Vandenburg,  and  Reno  are  in  good  agreement  with  the  isotachs.  The  two  charts  are 
in  very  good  agreement  as  to  the  general  location  of  the  height  contours  and  the 
closed  low  over  southern  California.  This  type  of  agreement  was  prevalent 
throughout  the  case  at  upper  levels. 

C.  LOWER  LEVELS  AT  1200  UTC  2  APRIL 

Examination  of  the  surface  charts  for  1200  UTC  2  April  yields  a  dramatic 
difference  compared  to  the  NCEP  chart.  Figure  23,  the  NCEP  chart,  has  a  broad 
inverted  trough  analyzed  that  extends  over  the  adjacent  offshore  region.  The 
COAMPS  surface  chart  (Fig.  36)  has  two  distinct  small-scale  troughs  analyzed,  one 
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in  the  lee  of  the  Sierras  and  one  over  the  coast  in  northern  California.  The 
automated  analysis  system  used  to  generate  NCEP's  surface  chart  does  smooth 
through  small-scale  features,  accounting  for  some  of  the  difference.  But  even  so,  the 
observations  are  not  dense  enough  to  confirm  the  small-scale  troughs.  Since  the 
troughs  are  associated  with  terrain  features,  reduction  error  might  be  playing  a  role. 
However,  this  SLP  analysis  utilizes  the  horizontal  reduction  techniques  proposed  by 
Mesinger  (1990)  which  tends  to  yield  less  reduction  error  than  conventional  vertical 
extrapolation  techniques  (Mesinger  and  Treadon  1995,  Pauley  1997).  It  is 
noteworthy  that  the  western  of  these  two  small-scale  troughs  is  roughly  below  the 
small-scale  500  mb  trough,  suggesting  that  this  feature  has  vertical  continuity. 

Further  analysis  of  the  COAMPS  surface  chart  reveals  the  interesting  features 
in  the  surface  isotachs  in  northern  California.  The  broad  windspeed  maximum  is 
located  off  shore,  similar  to  that  seen  in  Pauley  et.  al.  (1996).  Notice  again  the 
distinct  wind  maxima  that  occur  over  the  land  in  conjunction  with  the  small-scale 
troughs.  The  observational  network  is  not  sufficient  to  explicitly  verify  these  features, 
but  there  is  a  15  m/s  observation  in  north-central  California  in  proximity  to  the  wind 
maxima.  The  700  mb  COAMPS  chart  (Fig.  37)  shows  the  same  pattern  of  small-scale 
troughs  and  wind  maxima  as  the  surface  chart.  The  correlation  of  the  maxima  at 
these  two  levels  suggested  the  existence  of  some  type  of  disturbance;  the  location  of 
the  disturbance  in  the  lee  of  the  Sierra  suggests  a  mountain  wave. 
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VI.  GRAVITY/MOUNTAIN  WAVES 


A.  UPPER  FRONT  INTERACTION 

Figure  38  shows  the  500  mb  vertical  velocity  (w)  field  (shaded  at  5  cm/s 
increments)  and  the  height  contours  at  1200  UTC  2  April,  the  same  time  that  the 
three  separate  low-level  wind  maxima  are  the  most  dominant.  A  negative  value 
implies  downward  vertical  motion,  positive  upward.  There  is  clearly  strong 
downward  and  upward  motions  in  the  northern  California  region  downstream  of  the 
Sierras.  Upon  further  examination,  the  downward  vertical  motion  maxima  are  in 
approximately  the  same  location  as  the  lower  level  wind  maxima  that  were  discussed 
earlier.  Also  notice  the  vertical  motion  maxima  in  Baja,  California,  another  region 
that  has  a  complex  terrain  (Fig.  6).  This  vertical  motion  pattern  suggests  the  presence 
of  a  mountain  wave  as  it  would  be  depicted  in  plan  view.  Furthermore  the 
correlation  between  the  intense  westernmost  downward  motion  maximum  and  the 
small-scale  500  mb  trough,  suggests  that  the  trough ing  is  a  hydrostatic  response  to 
adiabatic  warming  associated  with  the  mountain  wave.  The  500  mb  temperature 
gradient  (Fig.  39)  agrees  with  this  premise.  Notice  the  strong  gradient  (shaded  at 
increments  of  3  K/lOO  km  starting  at  3  K/100  km)  associated  with  the  western  most 
downward  maximum.  The  isotherms,  represented  by  thin  solid  lines  at  2°C,  also 
show  tighter  packing  in  the  same  area.  This  in  turn  would  strengthen  the  jet  streak 
due  to  thermal  wind  considerations  thereby  enhancing  the  upper  front  as  well, 
representing  positive  feedback.  The  jet  streak  at  300  mb  does  in  fact  strengthen  from 
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0000  UTC  (Fig.  32)  to  1200  UTC  (Fig.  34)  2  April  as  can  be  seen  by  comparing  the 
two  charts.  The  increase  in  speed  would  also  tend  to  strengthen  the  secondary 
circulations  associated  with  the  jet  making  for  a  stronger  upper  front.  The  stronger 
the  upper  front,  the  stronger  the  stable  layer  which  implies  an  enhanced  mountain 
wave  completing  the  cycle. 

Figures  40  and  41  depict  a  cross-section  through  north-central  California  at 
1200  UTC  2  April  in  plan  view  and  the  vertical  view.  The  shaded  areas  are  the 
isotachs  at  5  m/s  increments,  the  dashed  lines  are  isentropes  and  the  dark  solid  lines 
represent  PV  Units  (PVU)  of  1.6  for  the  lower  and  3.0  for  the  upper  line.  The  1.6 
PVU  line  is  representative  of  the  tropopause  and  upper  front  signature.  Comparison 
of  this  cross-section  and  that  of  the  Boulder  windstorm  examined  by  Durran 
(1990)(Fig.  5)  show  remarkable  similarity.  Both  show  the  tight  packing  of  the 
isentropes  in  the  stable  layer  upstream  of  the  mountaintop  and  the  strong  wave 
pattern  in  the  isentropes  downstream.  The  tight  packing  of  the  isentropes  at 
tropopause  level  over  the  first  crest  of  the  mountain  wave  downstream  of  the 
mountain  peak  is  also  present  in  both  figures.  The  isotachs  in  Fig.  41  also  show  the 
wave  pattern.  This  field  is  not  depicted  in  the  Boulder  case  so  direct  comparison  is 
not  possible.  Oakland  and  Reno  rawinsonde  data  support  the  isotachs  but  cannot 
resolve  the  wavy  structure.  Notice  that  a  strong  upper  front  is  not  indicated  by  the 
PVU  contours  at  this  time  even  though  the  stable  layer  is  present  upstream  of  the 
mountain.  Additional  cross-section  analysis  is  necessary  to  resolve  the  influence  of 
the  upper  front  in  this  case. 
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A  cross-section  through  the  same  plane  but  at  0000  UTC  2  April  is  shown  in 
Fig.  42.  There  is  a  strong  upper  front  signature  defined  by  the  1.6  PVU  contour  that 
intersects  the  mountain  topography.  Notice  however  that  the  wave  pattern  is  not  as 
dominant  as  it  is  at  1200  UTC  2  April  (Fig.  41).  So,  at  this  time,  0000  UTC  2  April, 
the  upper  front,  a  stable  layer  by  definition,  interacts  with  the  Sierra  Nevada  range 
and  excites  a  mountain  wave  as  defined  by  Durran  (1990).  As  the  upper  front 
continues  to  propagate  with  the  jet  streak  from  0000  UTC  to  1200  UTC  2  April  (see 
synoptic  overview),  the  strong  frontal  signature  is  no  longer  present  in  this  cross- 
section  plane,  but  the  resultant  wave  remains  and  is  supported  by  a  stable  layer 
upstream  that  is  not  connected  to  the  dominant  upper  front. 
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VII.  CONCLUSIONS 


This  thesis  presented  a  case  study  of  an  upper  front  and  mountain  wave  event 
that  contributed  to  high  winds  in  central  California  on  2  April  1997,  This  chapter 
identifies  the  primary  results  of  the  study  derived  from  the  COAMPS  analyses. 

A.  UPPER  FRONT 

The  upper  front  associated  with  the  jet  streak  that  moved  through  northern 
and  central  California  on  2-3  April,  1997  excited  a  mountain  wave  just  north  of  Lake 
Tahoe,  as  the  front  interacted  with  the  Sierra  Nevada  range  in  that  area.  The  stable 
layer,  a  signature  of  an  upper  front,  was  near  mountain-top  level  and  on  the  upwind 
side,  a  necessary  situation  to  generate  mountain  waves.  The  wave  initiated  in  this 
region  extended  to  offshore  regions  west  of  San  Francisco.  The  wave  interacted  with 
the  surface  in  two  places  as  evidenced  by  the  two  distinct  wind  maxima  at  the 
surface  that  were  coincident  with  the  wave  pattern.  Thus  the  high  surface  winds 
generated  in  those  areas  were  at  least  in  part  due  to  higher  velocity  winds  aloft  being 
brought  to  the  surface  by  the  wave. 

Although  not  directly  associated  with  the  upper  front,  downslope  winds  were 
also  modeled  along  the  lee  of  the  Sierras  far  to  the  south.  The  impact  of  the 
downslope  winds  on  the  total  surface  wind  field  is  not  exactly  known.  Adiabatic 
warming  of  the  descending  air  would  reduce  the  surface  pressures  on  the  lee  side 
and  would  result  in  stronger  winds  due  to  the  tightening  of  the  pressure  gradient. 
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B.  DIFFERENCES  FROM  1991  "I  5"  DUST  STORM 

The  biggest  difference  between  this  case  and  the  1991  case  was  the  direction 
of  the  upper  level  winds  and  accompanying  jet  streak.  In  the  1991  case,  the  upper- 
level  winds  were  from  the  north-northwest  whereas  in  this  case  the  winds  and  jet 
streak  are  from  the  northeast.  The  northeasterly  flow  caused  the  jet  streak  and  upper 
winds  to  interact  with  the  Sierra  Nevada  range  producing  mountain  waves  and 
disrupting  the  flow.  The  1991  winds  were  nearly  parallel  to  this  mountain  range  and 
had  no  obvious  interaction  with  the  topography. 

A  second  difference  that  may  also  account  for  the  absence  of  mixing  upper- 
level  winds  to  the  surface  was  the  stability  of  the  lower  levels.  This  case  had  a  more 
stable  lower  layer  than  the  1991  case  which  may  have  inhibited  mixing. 

The  time  of  year  and  climatological  conditions  are  different  as  well.  The  1991 
dust  storm  occurred  in  early  November  after  a  period  of  unprecedented  drought  in 
California.  This  case  occurred  in  early  April  following  an  abnormally  wet  winter. 
The  degree  of  soil  heating  based  on  moisture  and  vegetation  would  be  less  for  this 
case  resulting  in  less  heating  of  the  boundary,  also  thereby  increasing  the  stability  of 
the  lower  levels  as  observed. 

C.  COAMPS  AND  MOUNTAIN  WAVES 

Because  COAMPS  has  only  been  in  use  for  a  short  time  and  is  not  yet 
operational,  there  is  not  a  lot  of  information  as  to  how  accurately  the  model  depicts 
meteorological  phenomena.  However,  one  case  study  was  done  that  analyzed  the 
ability  of  COAMPS  to  model  mountain  waves.  The  case  study  by  Doyle  et  al.  (1995) 
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describes  how  COAMPS  successfully  modeled  a  downslope  wind  and  mountain 
wave  event  that  occurred  in  Oppdal,  Norway  in  January  1995.  Although  in  this  case 
the  stable  layer  was  a  surface  warm  front,  cross-section  comparison  between  that 
event  and  the  April  1997  event  are  very  similar.  Figure  43  shows  four  panels 
depicting  the  Oppdal  wind  storm.  In  Fig.  43a,  isentropes  are  presented  with  the 
underlying  topography.  Compare  this  to  the  isentropic  pattern  for  the  April  1997 
case  (Fig.  41)  and  notice  the  obvious  similarity  in  the  theta  surface  pattern.  The 
Oppdal  cross  section  was  run  at  3  km  resolution  vice  27  km  which  accounts  for  the 
smaller  amplitude  waves  shown.  The  other  three  panels  of  Fig.  43  show  the 
geographic,  topographic  and  mesh  boxes  (Fig.  43b)  used  for  Oppdal,  the  cross 
section  of  horizontal  wind  speed  (Fig.  43c),  and  the  cross  section  plane  in  the  plan 
view  (Fig.  43d). 

There  are  differences  between  the  two  cases,  namely  the  vertical  position  of 
the  stable  layer  and  the  finer  mesh  used  in  the  Oppdal  case.  However,  the 
underlying  premise  that  COAMPS  has  been  known  to  successfully  model  mountain 
wave  events  cannot  be  refuted.  This  gives  one  greater  confidence  in  the  model 
output  for  the  April  1 997  case,  even  though  observations  remain  insufficient  to  verify 
the  event. 
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VIII.  RECOMMENDATIONS 


The  following  recommendations  are  provided  for  possible  follow-on  work  to 
this  case  or  similar  ones.  First,  seek  additional  non-conventional  observations  that 
may  be  available  for  the  area  of  concern  to  verify  the  model  output  for  this  case, 
especially  the  high  surface  winds.  Although  several  of  these  sources  were  examined, 
more  may  be  available  through  other  networks  or  organizations. 

Second,  run  this  case  with  a  higher  resolution  for  COAMPS.  At  a  9  km  mesh 
for  instance,  smaller  amplitude  features  may  be  analyzed.  In  addition,  a  100  m 
terrain  field  could  be  used  instead  of  the  20  km  used  in  this  study  due  to  the  finer 
resolution  of  the  inner  most  grid. 

Third,  compare  the  COAMPS  output  with  those  of  other  mesoscale  models 
such  as  the  10  km  eta  model  from  NCEP  or  the  MM-5,  run  locally  at  the  Naval 
Postgraduate  School.  Not  only  could  this  help  to  verify  the  COAMPS  output,  but  it 
would  also  provide  direct  comparison  of  a  non-hydrostatic  model  (COAMPS  or  MM- 
5)  to  a  hydrostatic  model  (10  km  eta)  and  their  ability  to  model  mountain  waves  and 
downslope  winds. 

Fourth,  run  COAMPS  forecast  fields  out  through  tau-24  or  tau-36  to  determine 
how  well  the  model  can  forecast  this  case  and  also  see  how  far  in  advance  the  model 
can  successfully  predict  this  event. 

Finally,  use  the  COAMPS  or  other  mesoscale  models  to  examine  other 
northeasterly  wind  events  over  the  Sierra  Nevada  range  and  see  if  they  also  yield 
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similar  patterns  to  what  occurred  in  this  case.  Differences,  if  there  are  any,  would 
be  useful  to  forecasters  to  better  their  prediction  knowledge  and  skills  for  mountain 
wave  and  downslope  wind  events. 


APPENDIX.  FIGURES 


HORIZONTAL  SHEAR 


Fig.  1.  Schematic  depictions  on  a  constant  pressure  surface,  of  straight  frontal  zones 
characterized  by  (a)  confluence  and  diffluence  associated  with  a  jet  maximum,  and 
(b)  horizontal  shear  in  the  presence  of  a  positive  along-front  thermal  gradient. 
Geopotential  height  contours,  thick  solid  lines;  isotachs  of  the  along-front  geostrophic 
wind  component  U,  thick  dashed  lines;  isentropes,  thin  dashed  lines.  Heavy-arrows 
indicate  the  jet  axis  in  (a)  and  the  sense  of  cross-front  shear  of  the  along-front  wind 
component  in  (b).  (from  Shapiro  and  Keyser  1990,  Fig.  10.6) 
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Fig.  2.  Schematic  illustration  on  a  constant  pressure  surface  of  various  idealized 
configurations  of  potential  temperature  and  along-front  geostrophic  wind  for  straight 
upper-tropospheric  jet  maximum.  Thick  solid  lines,  geopotential  height  contours;  thick 
dashed  lines,  isotachs  of  the  along  front  wind  component;  thin  solid  lines,  isentropes  or 
isotherms;  thick  solid  arrows,  sense  of  cross-front  ageostrophic  wind  componeriT at  level 
of  maximum  wind;  plus  and  minus  signs,  sense  of  mid-tropospheric  pressure-coordinate 
vertical  velocity,  (a)  Pure  confluence  and  diffluence  in  the  absence  of  along-je,tJhermal 
advection;  (b)  pure  horizontal  shear  with  80/ax >0  in  the  absence  of  the  effect  of 
confluence  and  diffluence.  (c)-(f)  Mixed  cases  of  confluence/d iffluence  and  horizontal 
shear:  (c)  along-jet  cold  advection,  (d)  along-jet  warm  advection,  (e)  jet  in  thermal  ridge, 
(0  jet  in  a  thermal  trough,  (from  Keyser  and  Shapiro  1986,  Fig.  23) 
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Fig.  3.  Schematic  diagram  of  the  present-day  model  for  analyzing  upper-level  fronts  and 

denote  frontal  and  tropopause  boundaries-  dashed  line  deoicts 
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Fig.  4.  Cross  section  from  the  1991  dust  storm  case  (Pauley  et  al.  1996)  Pressure 
IS  the  y-axis  and  distance  along  the  cross  section  is  on  the  x-axis.  Isentropes  are 

dashed  lines,  isotachs  are  shaded  at  5  m/s  increments,  and  the  bold  solid  lines  are 
1.6  (lower)  and  3.0  (upper)  PVU. 
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Tei^rain 


Fig.  6.  The  terrain  field  at  the  inner  27  km  resolution  domain  used  for  the  COAMPS 
model  run,  shaded  at  500  m  increments  starting  at  500  m. 
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'“'  °°°°  1  April  1997.  Height  contours  are 

solid  lines  (6  dm  increment)  and  isotherms  are  dashed  (5°  C  increment). 
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are  clashed  (5°  C  increment). 


Fig.  10.  NCEP  difax  850  mb  chart  for  0000  UTC  1  April  1997.  Height  contours  are 
solid  lines  (3  dm  increment)  and  isotherms  are  dashed  (5°  C  increment). 
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Fig.  12.  Same  as  Fig.  7  but  at  1200  UTC  1  April  1997, 


Fig^  13.  Same  as  Fig.  8  but  at  1200  UTC  1  April  1997. 
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Fig.  15.  Same  as  Fig.  11  but  at  1200  UTC  1  April  1997 


Fig.  16.  Same  as  Fig.  7  but  at  0000  UTC  2  April  1997. 
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Fig.  17.  Same  as  Fig.  8  but  at  0000  UTC  2  April  1997. 


Fig.  18.  Same  as  Fig.  10  but  at  0000  UTC  2  April  1997, 
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Fig.  19.  Same  as  Fig.  11  but  at  0000  UTC  2  April  1997 


Fig.  20.  Same  as  Fig.  7  but  at  1200  UTC  2  April  1997. 


i*. 


I  ’  .  •Vv.vV^iiU  >rriv1 

fMl!fA;:i4ii^Kll&fediM^ 


lii  lL'^i  .! 


h»  OJini 

iFf  H  I  ■  \  i-t::  'A  ■  /oX 

niiili 

h-iftfUwii—ifUUi ,IK  i .• '-Sts tev 


^Ki  nNini 


\  ;  iss  'mihX-^r^ 
\  .’  r  i-'W  >  u» 

=T/  L  .  ■  .  i  ••■<4“  ii'iaii  .ys*  >•- 


Fig.  *21.  Same  as  Fig.  8  but  at  1200  UTC  2  April  1997. 


Fig.  22.  Same  as  Fig.  10  but  at  1200  UTC  2  April  1997, 
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Fig.  24.  Same  as  Fig.  7  but  at  0000  UTC  3  April  1997. 


Pig.  25.  Same  as  Fig.  8  but  at  0000  UTC  3  April  1997. 


Fig.  26.  Same  as  Fig.  10  but  at  0000  UTC  3  April  1997. 
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Fig.  27.  Same  as  Fig.  1 1  but  at  0000  UTC  3  April  1997. 


Fig.  28.  COAMPS  300  mb  chart  at  0000  UTC  1  April  1997.  Height  contours  are 
solid  lines  (12  dm  increment),  isotachs  are  shaded  at  5  m/s  increments  starting  at  30 
m/s.  Observations  are  in  m/s:  flag  =  25  m/s,  full  barb  =  5  m/s,  and  half  barb  =  2.5 
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Fig.  29.  COAMPS  surface  chart  at  0000  UTC  1  April  1997.  Isobars  are  soRd  lines 
(4  mb  increment),  isotachs  of  the  average  surface  wind  speed  are  shaded  in  5  m/s 
increments  starting  at  15  m/s.  Observations  are  in  m/s:  full  barb  =  5  m/s  and  half 
barb  =  2.5  m/s.  The  average  windspeed  here  is  defined  as  the  average  of  the  lowest 
four  sigma  levels,  corresponding  to  the  lowest  90  m.  Values  are  weighted  by  the 
layer  depth  in  forming  the  average. 


QOZ  2  Apr  97  300  mb  windspeed 


Fig.  32.  Same  as  Fig.  28  except  for  0000  UTC  2  April  1997 


Fig.  35.  COAMPS  500  mb  chart  at  1200  UTC.2  April  1997.  Height  contours  are 
solid  lines  (6  dm  increment).  Other  fields  as  described  in  Fig.  28. 
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12Z  2  Apr  '97 _ 700  mb  windspeed 


Fig.  37.  COAMPS  700  mb  chart  at  1200  UTC  2  April  1997.  Height  contours  are 
solid  lines  (3  dm  increment).  Other  fields  as  described  in  Fig.  28. 
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12Z  2  Apr  ^97  500  mb  model  vertical  velocity 


Fig.  38.  COAMPS  500  mb  chart  at  1200  UTC  2  April  1997.  Height  contours  are 
solid  lines  (6  dm  increment).  Vertical  velocity  field  (w)  is  shaded  in  10  cm/s 
increments  beginning  at  5  cm/s.  Negative  values  indicate  downward  motion  and 
positive  values  indicate  upward  motion. 
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Fig.  39.  COAMPS  500  mb  chart  ai 
solid  lines  (6  dm  increment).  Tem| 
Temperature  gradient  is  shaded  in 


Rg.  40.  COAMPS  surface  chart  at  1200  UTC  2  April  1997.  Isobars  are  solid-lines 
(4  rob  increment).  Cross  section  location  denoted  by  line  Z  to  Z'.  10  m  isotachs  are 
shaded  in  5  m/s  increments  starting  at  10  m/s. 
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Fig.  41.  COAMPS  cross  section  at  1200  UTC  2  April  97  through  the  Z  to  Z' 
Identified  in  Fig.  40.  Cross  section  details  same  as  in  Fig.  4. 
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Fig.  42.  Same  as  Fig.  41  except  for  0000  UTC  2  April  1997. 
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Fig.  43.  Diagrams  from  Doyle  et  al.  1995.  a)  is  isentropic  cross  section,  b)  is 
topographic  and  COAMPS  grid  map,  c)  is  cross  section  of  horizontal  wind  speed  in 
m/s;  d)  is  plan  view  of  cross  sectional  plane  C  to  C'. 
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